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Abstract 

Several proposals for quantum computation utilize a lattice type architecture 
with qubits trapped by a periodic potential. For systems undergoing many body 
interactions described by the Bose-Hubbard Hamiltonian, the ground state of the 
system carries number fluctuations that scale with the number of qubits. This pro- 
cess degrades the initialization of the quantum computer register and can introduce 
errors during error correction. In an earlier manuscript we proposed a solution to 
this problem tailored to the loading of cold atoms into an optical lattice via the 
Mott Insulator phase transition. It was shown that by adding an inhomogeneity to 
the lattice and performing a continuous measurement, the unit filled state suitable 
for a quantum computer register can be maintained. Here, we give a more rigorous 
derivation of the register fidelity in homogeneous and inhomogeneous lattices and 
provide evidence that the protocol is effective in the finite temperature regime. 

1 Introduction 

The goal of building a quantum computer has spurred tremendous progress in coher- 
ent control and measurement of small quantum systems. In order to fully realize the 
promised computational speedup of a quantum device, the underlying system should 
be scalable to a large number of information carriers or qubits. Indeed, the first two 
criteria delineated by DiVincenzo [1] for scalable quantum computation are: 

• A scalable physical system with well characterized qubits 

• The ability to initialize the state of the qubits to a simple fiducial state 

In many systems, the first criterion can be met by increasing the number of storage 
components for the qubits, e.g. in solid state systems the number of dopant qubits 
in the bulk material could be increased, or for ions large scale micro-trap arrays have 
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been proposed [2]. There are two main approaches to satisfy the second criterion [1]. 
One is to allow the system to interact with the environment and "naturally" cool to 
its ground state and thereafter use this state as the initial state. The other is to actively 
cool each qubit by projective measurement to a fiducial state |0). A problem arises with 
these approaches when the ground state of the many body Hamiltonian is not a suitable 
initial state. This is the case for bosonic qubits embedded in systems with periodic 
confinement, for example in josephson junction arrays [3] and neutral atoms trapped in 
electromagnetic microtraps [4] or optical lattices [5] . For these systems, the underly- 
ing dynamics is Bose-Hubbard like and the ground state contains residual coherences 
described by non-zero number fluctuations in each mode. We show in Sec. |2]that these 
fluctuations scale with the number of qubits so that if not corrected, the dynamics will 
impart a constraint to scalability. In many cases, the tunneling energy between lattice 
sites is weak so that the overall error probability per qubit is small, however, if the 
system is projected onto a local number state basis during error correction cycles, the 
error can become appreciable [6]. 

We propose a solution to this problem which reduces number fluctuations to neg- 
ligible size even for large numbers of qubits. The technique is exemplified in a three 
dimensional optical lattice filled with cold neutral bosonic atoms. Optical lattices are 
three dimensional potentials created by counter-propagating beams of laser light that 
trap atoms by the conservative electric dipole potential [7]. This is a particularly attrac- 
tive system to study because it has been theoretically shown [8], and recently demon- 
strated in the laboratory [5], that an optical lattice can be loaded from an atomic BEC. 
If the lattice is turned on adiabatically with respect to low lying many body excitations, 
then the superfluid-like BEC undergoes a phase transition to the Mott insulator (MI) 
state characterized by predominantly the same number of atoms in each lattice well [9] . 
A difficulty with using this mechanism to prepare a unit filled register of atomic qubits 
arises due to imperfect filling of the lattice as well as number fluctuations intrinsic to 
the Bose-Hubbard dynamics. 

In [10] we described how to eliminate number fluctuations in the many body ground 
state. The approach was to first introduce a quadratic magnetic trap that acts to fill 
sites with missing qubits in the center, and second to project out components of the 
many-body wavefunction with multiply occupied lattice sites. In this paper we give a 
more rigorous derivation of the dynamics due to first order couplings outside the unit 
filled state in both homogeneous and inhomogeneous systems. The essential result is 
that deep in the Mott Insulator regime, the dynamics can be effectively modeled with 
a restricted basis involving two level couplings between the unit filled state and first 
excited states with double occupancy in one lattice site. We then show in Sec. |3lthat by 
turning on an external field components of the many body state with multiply occupied 
wells can be projected out. Provided the measurement strength is sufficiently large, the 
system does not evolve out of the restricted basis and the measurement can maintain a 
unit filled register for the duration of quantum computation. Finally, we give evidence 
that the proposal can be made robust even for finite temperature systems. 
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2 Bose-Hubbard Dynamics 



The Bose-Hubbard Hamiltonian describes a system of bosons on a lattice having 
a kinetic energy associated with tunneling between nearest neighbor lattice sites and a 
pairwise interaction at each lattice site: 

hbh = -Y. +Le(j)«j + - 1) (1) 

<ij> j 

Here aj are the bosonic annihilation operators and nj — aj aj the number operators for 
an atom in the lowest vibrational state of lattice well j = {ji,j2....,jd) with d the dimen- 
sionality of the lattice. The energy offset at each lattice site is £(j) which models a con- 
tinuously varying external potential. Henceforth we assume that the energies J and U 
describing tunneling and on-site interaction energies respectively are site-independent. 
The notation < i,j > indicates that the sum is over nearest neighbors. In the tight- 
binding model, the nearest neighbor tunneling energy J is defined as one fourth the 
band width of the lowest occupied band. The dynamics of an ultracold bosonic gas in 
an optical lattice loaded in such a way that only the lowest vibrational state is occupied 
can be described by Hbh where the system parameters are controlled by laser light. 
In a real optical lattice, tunneling in one, two, or three dimensions can be achieved by 
tuning the laser intensities and detunings so that dynamics can be effectively frozen in 
some directions. We assume a square lattice with tunneling dynamics in d dimensions 
through a separable potential barrier V{x) = VY4^iCos^{kxi). The tunneling rate is 

closely approximated by J/h = A/ {s/Tih)ER{V /ERfl'^e'^^/^^ [11], where the recoil 
energy is Er — {hk)^/2m (m =atomic mass). The on site interaction is a result of the 
ground state collisions described by the i-wave scattering length Oj between two atoms 
each in the motional state (|)(x) and is given hy U — / iix|(|)(x) |^. 

It was recognized early on that the MI transition might be an efficient way to ini- 
tialize a register of atomic qubits in an optical lattice for use in quantum information 
processing. A key advantage of loading from a BEC is the availability of an initially 
high phase space density which can be frozen to the MI state with atoms occupying 
most lattice sites. For the homogeneous system (e(j) = 0), only commensurate fillings 
give rise to a MI transition. For the purposes of quantum computation, one particle 
per well is desirable. In practice, this is difficult to achieve directly because the pre- 
cise number of atoms is unknown and the lattice strength is not perfectly uniform on 
the boundaries. There are proposals to prepare unit filled lattices using dissipative 
techniques involving filling the lattice atom by atom [12] or using Raman side-band 
cooling [13]. Additionally, it has been shown that one can repair imperfect filling from 
a BEC via an adiabatic transfer mechanism between two sublevels of each atom [14]. 
While these techniques can efficiently initialize the lattice, the unit filled register is not 
a stationary state of the system and the dynamics resulting from the inter-well coupling 
degrades the register fidelity defined as the population in the unit filled state. In Sec. 
I2.1l we discuss the dynamics in a homogeneous system assuming that an initial state 
with unit filling has been prepared. We show that in the strong coupling limit, J/U ^ 1, 
the dynamics can be restricted to a subspace whose dimension grows linearly with the 
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number of particles in the system. The fidelity, in this limit, is then derived explicitly. 
In Sec. I2.3l we show that an external quadratic potential can be used to help prepare a 
unit filled state in an appropriate subspace of the lattice, which will constitute a quan- 
tum computer register. The dynamical fidelity derived in the homogeneous case can 
then be mapped to the one in the register. 



2.1 Homogeneous approximation 

In order to understand the relevant dynamics in our system we first discuss the prop- 
erties of an idealized homogeneous lattice with periodic boundaries and unit filling, 
N = M where is the number of atoms and M is the number of wells. The regime of 
interest is the strong coupling limit, which is experimentally achievable in an optical 
lattice because the tunneling decreases exponentially with the trap depth. To zero or- 
der, neglecting the kinetic energy, the ground state of the system is the unit filled state 
\T) = rij flj |0) which is the Fock state with exactly one atom per well and energy zero. 
A small but non-zero J /U introduces number fluctuations in the ground state which 
provide a small residual coherence across the system [15, 16]. We can calculate the ap- 
proximate ground state by including states which couple to \T) to second order in Hbh- 
The first A^(A^ — 1 ) degenerate excited states are Fock states with two particles in one 
well, one hole in another well and unit filling in every other site. We define such states 
particle-hole excitations and their energy is equal to U . To first order, |r) couples only 
to a state with a doubly occupied well neighboring the hole. Because Hbh is invariant 
under translations of the lattice, this state, when properly symmetrized, is 

\si)= y ■ (2) 

To second order, the coupling is to three symmetrized states: the state with two particle- 
hole excitations |Cii), the state with three particles in one site and holes in two neigh- 
boring sites \Qi), and 1^2), the symmetrized state where the doubly occupied site and 
the hole are next nearest neighbors. The former two states are given bj 

=/(A?,£/)£<i j>_i<pq>fl?'fljaqflp|r), \Qi) =g(A?,£/)i:<ijk>fliflj- flk|r), and have 
energies 2U and 3f/ respectively. The second order coupling is proportional to the 
inverse of the normalization factors of the three states: for |Cii), the normalization 
is f[N4) = {2dN{dN -3-A{d-\)))-'^l^-, for \Qi) the normalization is g(N,d) = 
{6d{2d- \)N)-^I^, and for \S2) it is {AdN)-^l^. Provided dN < {U/jf, then the per- 
turbation expansion is convergent and to first order the ground state is approximately, 

I'Pg) = a[\T)+2^-VdN\Si)^ (3) 

The normalization constant is a = (1 +4Nd{J /U)^)^^^^ . To second order, the ground 
state energy is £ = —AdNJ^/U. Notice the fidelity F ~ = decreases 

linearly with the number of atoms for fixed J /U . 
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2.2 Homogenous dynamics 



Any mechanism used to prepare a register of qubits in the unit filled state will suffer a 
degradation in fidelity, because the state \T) is not a stationary state of the system. In 
this section we derive the time dependent fidelity in a homogeneous commensurately 
filled system prepared at time f = in the state | T) . While our results can be general- 
ized to higher dimensions, in the remainder of this section we assume tunneling along 
one dimension only. To compute the first order coupling outside of \T) we need to 
rediagonalize Hbh in the subspace spanned by the states with one particle-hole pair. In 
principle we should consider all N{N — 1) states. However by translational invariance 
of Hbh, the dynamics is restricted to the lN/2\ invariant states in the subspace, where 
[jcj denotes the greatest integer less than or equal to x. The states in the particle-hole 
subspace are represented as 



[N/2\ 

n=l l'-il=" 



(4) 



In order to diagonalize the Hamiltonian the coefficients .sjj must satisfy the following 
recurrence relation 



- 37«,+i + 1 ) = (£. - U)/,- si = 0, = /„ . 
The solutions of the above equations can be shown to be given by 

'%{2r- 1) 



2 



7r(2r- 



N 



Ey = U — 6/ COS 
The state vector at any time can be written as 

mt))^c,{t)\T 



1,- 



1, 



XN/2\ 



yNi2\ 



with co(0) = 1. 



(5) 

(6) 
(7) 

(8) 



Solving the Schrodinger equation to first order in perturbation theory, the time de- 
pendent fidelity to be in the unit filled state F{t) — |co(f)P can be estimated to be 



F{t) 



[NI2\ 

1- E MOI' 

r=0 
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J 
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E, 



7t(2r+i: 



sm 



E,-t/h 



(9) 



(10) 



The sum in Eq. [Ojcan be solved in the strong coupling regime of interest and the closed 
form expression is 



F{t) ==1-8 



1 



cos{ut/h)Ji {6tj/n) 
3Jt/h 



(11) 
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where Jn is the nth Bessel function of first kind. This solution compares well with 
exact numerical simulations for small number of wells [N <%) and for larger in the 
subspace with at most one site with two atoms. The behavior of F{t) is determined 
by fast oscillations with frequency equal to one, modulated by longer oscillations with 
frequency determined by the zeros of Jy . For short times, Jt jfi^ 1 , the fidelity is 

^(0-l-16(^) W(^), (12) 

which corresponds to the Rabi oscillations of an effective two level system spanned by 
|r) and \S\) . For later times the coupling to the other states become important. Notice 
that fidelity when time averaged over many oscillation periods Inh/U is {F{t)) = 1 — 
S{J/U)^N. Consequently, the deviation from a perfectly prepared register, described 
by 1 — {F{t)), is twice as bad as if the system were prepared in the stationary ground 
state \Wg) of the Bose-Hubbard Hamiltonian. This indicates that a lattice filled with a 
commensurate number of atoms in the Mott-Insulator state may create a more robust 
quantum computer register than one prepared dissipatively. 

2.3 Dynamics in presence of the external trap 

As mentioned in Sec. |2]in practice it is difficult to prepare an optical lattice with 
a commensurate number of atoms and wells. It can be arranged so that < M in 
which case we propose to use an inhomogeneous lattice with open boundaries cre- 
ated by a weak quadratic magnetic trap. The addition of the trap acts to collect atoms 
near the potential minimum and leaves empty wells (holes) at the edges. In exper- 
iments where the optical lattice is loaded from a EEC, the external trap is already 
present to confine the condensate. For simplicity we assume a spherical trap with 
oscillation frequency in each dimension given by coj. The magnetic confinement in- 
troduces a characteristic energy scale 5 = m/2{%/k)^(S]ij so that e(j) = 5^f^i jj. In 
order to insure that the on site interaction energy U is larger than the trapping en- 
ergy of the most externally trapped atoms, we require that the trap strength satisfy 
U > 5{d/2)\N^^'' /%''/^. Multiple atom occupation is therefore inhibited in any well 
in the ground state configuration. The register is defined by the subspace compris- 
ing K < N wells in the center region of the the trap (see Fig. [Q. The barrier space 
flanking ^ will act to suppress percolation of holes from the edges to the center. The 
probability for holes in due to tunneling through the barrier is estimated by comput- 
ing the product of transition probabilities arising from first order couplings from the 
boundary of the register at radius r„„„ to the outside edge of the occupied lattice at ra- 
dius r,nax- Accounting for tunneling from any site at the surface of the hyper-spherical 
shen at w, we find p,, ^ n[i|=V„,,/'1?,«i"''/^-/^(5(2|j| + l))-y{d/2)l - [(r[r,„,„ + 
l]/r[r,„„., + 2])^dri„a.^%''/^/{d/2)\] [7/(25)](2'-™'-2'-mm+2). The percolation probability 
Ph is negligible provided the barrier region is sufficiently large and J/ (2r,„/„5) < 1. 

Hereafter, we restrict our attention to the dynamics of the reduced state of the reg- 
ister: p(f ) = Tr^^j^ obtained by the tracing over spatial modes outside the register 
subspace of the entire many body wavefunction a. As with the analysis of dynamics 
in the homogeneous case, for simplicity we consider tunneling in one dimension only. 
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Figure 1 : Schematic of an inhomogeneous lattice filled with qubits with an onsite 
interaction energy U . An externally applied trapping potential of strength e(y) = 87'^, 
e.g. due to a magnetic field, acts to fill gaps in the central region of the trap. The center 
subspace of the lattice defines the quantum computer register containing K < N 
qubits. 



The degree of inhomogeneity is quantified by the ratio 5/7. For < 5/7 < 1, the en- 
ergy splitting between Fock states describing particle-hole pairs is small and the model 
of homogeneous dynamics is valid. For 5/7 > 1, to first order in perturbation theory, 
the state space is spanned by the unit filled state in the register subspace, \T),^, and the 
2K nearest neighbor particle-hole pairs \Sj): 

1— r + I did +1 d -i^a 

in^- n a]\0), \Sp.^ = ^^\T)^, \SJ)^^^\T)^. (13) 

j=-{K-l)/2 

Here we have introduced coordinates with the site j = coincident with the trap min- 
imum. For each j the states iSj")^ are distinguished by the two energetically dis- 
tinct orientations of a doubly occupied site and its neighboring hole with energies, 
= j/(lipA(2j-l)). Wedefinethe zero of energy coincident with the state | T) . 
We proceed to estimate the fidelity to be in the unit filled state in the register: 
■fk, (0 = Tr[a(f)|r)3(^3(^ (r|]. With an eye to the projective measurement, it is important 
to understand the free dynamics when the trap is present and the register is initialized in 
the pure state, p(0) = \T)^^ {T\. The quantity (f ) is generally difficult to compute 
because atoms couple into and out of the register. We can solve the simpler problem of 
the fidelity to be in the target state of a commensurately filled inhomogeneous lattice. 
This corresponds to a situation where the register occupies the full extent of a lattice 
{K = N) with periodic boundary conditions. In this case, to first order in J/U, the time 
dependent state can be written: 

\^>it))^CTit)\T)^+Y,e-'^^'P/'c,.{t)\Sf)^. (14) 

The fidelity Fca,„{8,K,t) = \ {T\W{t)) p can then be solved as in the homogeneous case. 
For 5 <C t/ and assuming the A^-body state is prepared in \T) it is given by 

Fcom{N,t) = l-i{J/uf{N -co&{Ut/Ti){l+&m{b{N -l)t/h)/iim{bt/h))). (15) 
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This solution compares well with exact numerical simulations for a small number of 
atoms (A^ < 8) and simulations for larger commensurate lattices with infinitely high 
boundaries using a restricted basis set of dimension A^(A^ — 1) + 1 consisting of the 
target state and all particle-hole pairs. Notice that for 5f — > 0, we recover the short time 
fidelity in the homogeneous case, Eq. [2] In the incommensurate case we can bound 
the fidelity F,j^ (f ) inside the register. Provided N > K the following inequalities on the 
time averaged fidelities hold: 

(F,„,,,(A^,f|a(0) = |r)(r|)) < (F^(f|p(0) = |r)^^(r|)) < (F,„,„(^^ 

(16) 

The lower bound arises because the probability to be in the unit filled state of a large 
commensurately filled lattice, given an initial state which is unit filled, is always less 
than or equal to the probability to be unit filled over a smaller subspace of K < N of 
a non commensurately filled lattice whose register is prepared in the unit filled state. 
This inequality holds provided the probability for holes to tunnel into the register ^ is 
small over relevant time scales. The upper bound is a consequence of the fact that unit 
filling in the register state is degraded because particles can tunnel in and out of the 
register. Therefore its fidelity is less than that of a commensurately filled lattice of the 
same size K prepared in the unit filled state. 

3 Measurement 

In Sec. I2.3l we showed that the system dynamics can be restricted to a set of two level 
couplings {\T)^ — > I'S'j^)}- We now sketch how to perform a continuous measurement 
to drive the register into the unit filled "target" state \ T)^. The full details are contained 
in [10]. The idea is to apply an external control field that is resonant with a coupling 
between the "faulty" states {I'Sj^)} and a set of excited states {|Mj^)j^}. If population in 
the excited states is easily measurable, for instance by the emission of photons during 
decay, then the presence of population in the particle-hole states can be monitored. It 
is vital that the coupling field be able to spectroscopically resolve the measurement 
transition without exciting the target state. This is possible if the atoms in multiply 
occupied wells see a shifted excited state, E{M^) —E{S'^) +E — U, where for example 
E is the energy of a dipole-dipole molecular state. The free energy of the system is 

Ho+Hbh = Lj.±E{Sf)\Sf)^^{Sf\+E{Mf)\Mf)^^{Mf\ 

~V2JLM\Sf)^^{T\ + \T)^^{Sf\). "-"^ 

When the control Hamiltonian with Rabi frequency £2 is turned on resonant with the 
measurement transition, the interaction Hamiltonian in the rotating frame is: 

Hi = LjAi\Vc\+EiSf))\Sf)^^{Sf\) + m]+EiSf)-U)\Mf)^^{Mf\ 

-V2Ji\Sf)^^{T\ + \T)^^{Sf\)+h£l/2i\Mf).j,.j,{Sf\ + \Sf)^^{Mf\)), 

(18) 

where the energy \Vc] includes possible off resonant energy shifts on atoms in singly 
occupied wells. Including decay from the excited states at a rate y, the non-unitary 
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Figure 2: Population in the unit filled register state \T).j^^ during continuous measure- 
ment of the register beginning in the Bose-Hubbard ground state Y^g). The plots show 
dynamics appropriate to tunneling in one dimension with U /J = 500. (a) Quantum 
trajectories corresponding to a null measurement result for three different register sizes 
K. The time scale to saturate the target state is independent of the number of qubits: 
tsat ~ K^^ (b) Long time dynamics for K — 5Q\,N = 55\ and finite detector efficien- 
cies r|. The population in \T)^ for r| = 1 is indistinguishable from one. Also shown is 
the oscillatory dynamics at fundamental frequency U described by Eq. E]if the mea- 
surement is turned off after the target state is reached. The arrow indicates p7- 7-(0). 



dynamics in the register is given by the non trace preserving master equation: 

p = -//r7[///,p]-Y/2^(|M±)^^(M±|p + p|M±)^^(M±|). (19) 

i.± 

Assuming low saturation of the excited states, the dynamics in the ground state is 

9s^T = -'Ps±r(^(^f) + l^c|)//i + /(pr,r-Ps±j±)V27/ri-p5±^K 

} J ■' J J J 

PTJ = iiPs±,T-pT.sf)^J/^ (20) 

Psf.sf = -'iPsf.T^PT.sf)^/^''~^Psf.Sf^- 

Here the non-selective measurement induces phase damping of the coherences at a 
rate K = H^y/ (8((f/ /Ti)^ + (y/2)^). When the environment is monitored, for example 
by watching for photons, the population in the particle-hole states will change. In 
the case of a null result, one obtains knowledge that the state is more likely in the 
target state and over time, in the "good" measurement regime [17], the probability 
to be in \T),j^ converges to one. The "good" measurement regime is given by the 
condition il/y ^ 1 < hK/2\/^J, which requires that the control field be strong enough 
to damp coherences on the time scale of the Bose-Hubbard dynamics but weak enough 
to not saturate the excited states. Simulations of successful register preparation and 
maintenance by measurement are shown in Fig.|2] 
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Figure 3: Equilibrium fidelity to be in the unit filled register state |r)j^ as a function 
of temperature in presence of the trap forM = \\,N ^9,K = 5. The relevant energies 
are = 60 and 5// = 3.375. Here w 6.0J/kB and T/, « IQJ/ks. The dashed line 
indicates the scaled energy Tciks/J where F{Tmux) ~ 0.33988. The fidelity at T = is 
F(r =0) w 0.99208. 



3.1 Measurement at finite temperature 

Up to now we have been focusing on dynamics of pure states of the many body sys- 
tem. The selective measurement is an entropy decreasing map, because it damps 
amplitude in multiply occupied wells. Therefore it can be effective even for mixed 
states at finite temperature. We hereby assume a thermal distribution of the eigen- 
states of the BH-Hamiltonian. The overall effect of finite temperature is to increase 
the weight of reduced Fock states other than the target state and therefore to decrease 
the fidelity. For a temperature Td ~ (t/ — e(^_i)/2)) /^b states are populated which 
have more than one particle in one or more sites in the register. The projective mea- 
surement is effective on these states and the primary consequence of their presence 
is to reduce the initial fidelity of the system. We can identify another temperature, 
7), « (£(A'+i)/2 ~ £((A'-i)/2)) /^B corresponding to the energy of the Fock state with a 
hole at site (K—l)/!. For T ^ T/, there is appreciable population in eigenstates which 
have holes in the register. The measurement is insensitive to population of states with 
holes in the register, and the probability that the measurement does not project the sys- 
tem into the target state, given a null measurement result, is then at least l/e when the 
temperature is greater than T/,. As a low initial fidelity can cause quantum jumps dur- 
ing the meaurement, large population of states with multiple occupancy in the central 
sites should also be avoided. This leads to the rough estimate T,„ax ~ min{Td, 7),}. The 
exponential growth of the Hilbert space with particle number makes the computation 
laborious. We have obtained numerical results for a model system consisting of 1 1 
sites and 9 particles by exact diagonalization of Hbh in a Hilbert space of dimension 
92378. The fidelity is shown in Fig.|3l The coupling ratio is chosen to he U /J = 60 
and the register is defined by the central 5 sites. In order to suppress tunneling of holes 
into the register, the ratio 5/7 has been chosen to be very large, 5/7 « 3.375. This does 
not correspond to a typical experimental situation, as it implies that the number of sites 
in the register for which e(y') <U is small (9, in this case). The fidelity drops to values 
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lower than F < 1/e at T^/ « 6.QJ /Ub — min{Tii,Th}, see caption of Fig.|3] For differ- 
ent setups Td and T/, can assume approximately the same value. In the experimentally 
relevant setup discussed in relation to the measurements, where = 551 and = 501, 
Td w 50J/kB, and T,, w mj/ke. 

4 Summary and conclusions 

In summary, we have proposed a protocol to eUminate number fluctuations in the many 
body ground state of a system undergoing Bose Hubbard dynamics. For concreteness 
we have focused on preparation of a register of cold atomic qubits in an optical lattice. 
By considering the free dynamics in both homogeneous and trapped systems we have 
shown that there is a loss of fidelity to occupy the unit filled register that scales with the 
number of qubits because this is not a stationary state. Because our analysis is based on 
first order perturbation theory we established the parameter regime where our analytic 
approach is valid by considering second order corrections. The dynamical degradation 
of the fidelity once the register is prepared is an issue that has to be considered not 
only in our continuous measurement proposal but also in other dissipative techniques 
proposed for unit filling initialization. We have suggested one approach to correct 
for this by performing a continuous measurement on the system. Numerical studies 
on small sized systems indicate that this protocol can be made robust even for finite 
temperature systems which will be important for real experimental implementations. 
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